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Abstract 

o . 

Q" 1 We investigate the effects of the accretion of phantom energy with non-zero bulk viscosity onto a 

Morris- Thorne wormhole. We have found that if the bulk viscosity is large then the mass of wormhole 
6J[), increases rapidly as compared to small or zero bulk viscosity. 
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Exotic tunnel like topology in spacetime commonly called wormhole arises as a solution 
to the Einstein field equations. A typical two mouth wormhole joins two arbitrary points 



either of the same spacetime or two different spacetimes. Morris and Thorne pQ suggested 
the existence of exotic matter for the stability of a wormhole. The stress energy tensor 
of the exotic matter must violate the Null energy condition {T^ u u^u u > 0) where 
m m is the future directed null vector. They concluded that an advanced civilization can 
produce a wormhole for interstellar travel by injecting sufficient amount of exotic matter 
in it. Recent interest in wormhole has arose due to the discovery of exotic phantom 
energy driving the accelerated expansion of the universe [21 E]- It has been proposed 
that wormhole can be stabilized by the accretion of phantom energy and it can result 
in increasing the size of the wormhole to engulf the observable universe [J]. We here 
consider a similar scenario where we incorporate the effects of bulk viscous stress in our 
calculations. The bulk viscosity is quite relevant in physical cosmology as it can cause 
expansion of the universe due to its negative pressure [5J. The presence of viscous fluid 
can also explain the observed high entropy per baryon ratio in the universe [6]. The 
formulation of the paper is adopted from [U [7] . 

The plan of this paper is as follows. In the second section, we present the relativistic 
model of accretion of viscous phantom energy onto a stationary wormhole. In third 
section, we consider two special cases of bulk viscosity in our model. Finally we conclude 
our paper. 
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2 Accretion onto wormhole 



We consider a stationary and spherically symmetric wormhole specified by the line ele- 
ment: 

ds 2 = -e*W + + r 2 (d6 2 + sin 2 6d<p 2 ). (1) 

The above metric has the following properties [H[9]: (1) The redshift function $(r) must 
be finite for all values of r thus no horizon exists outside the spacetime, (2) the shape 
function b(r) must satisfy the following conditions at the throat r = r Q : b(r ) = r Q and 



b'{r ) < 1, (3) b(r) < r for r > r and (4) the spacetime is asymptotically flat i.e. 







as r 



oo. 



Outside the wormhole, we assume the spacetime to be Friedmann-Robertson- Walker 
containing only one fluid, namely the phantom energy with non-vanishing bulk viscosity. 
The fluid is assumed to fall onto the WH horizon in the radial direction only which 
is in conformity with the spherical symmetry of the WH. Thus the fluid four velocity 
is = (w*(r), u r (r), 0, 0) which satisfies the normalization condition u^u^ = 1. The 
corresponding stress energy tensor for the exotic phantom energy is 

T^ = (p + p eff )uV+p eff g^. 



(2) 



Here p e ff = p+p V i SC , where p is the isotropic pressure and p visc = —3£,H is the bulk viscous 
pressure with £ is bulk viscosity. Using the energy momentum conservation T% u = 0, we 
get 

b(r 



-i 



ur 2 M- 2 (p + p eff ) 1 



b{r) 



- 1 = c l5 



(3) 



where u r = u = dr/ds is the radial component of the velocity four vector and C\ is a 
constant of integration. The second constant of motion is obtained by contracting the 
four velocity with the energy momentum conservation equation u^T^ v = 0, which gives 



w 



7V + (P+Peff, 



ur 



0. 



Integration of Eq. (4) gives the second constant of motion 



2 -2 



ur z M 



b(r) 



cxp 



ph 



dp' 



ft + Peff (ft) 



-A, 



(4) 



(5) 



where A is a constant of integration. Also ph and p^ is the energy density of the phantom 
fluid at the horizon of the WH, and at infinity respectively. From Eqs. (3) and (5) we 
have 



(P + P 



eff) 



1 - 



b(r) 



-1/2 



b{r) 



1 exp 



dp' 



P'+Peff(p') 



Co 



(6) 



here C 2 = —C\j A = A^p^ + p(p oa )). In order to calculate the rate of change of mass 
of WH M, we integrate the flux of the bulk viscous phantom fluid over the entire WH 
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horizon i.e. 

M = j>T r t dS. (7) 

Here T[ determines the energy momentum flux in the radial direction only and dS = 
y/—gd6dip is the infinitesimal surface element of the WH horizon. Using Eqs. (3) - (7), 
we get 

dM , / bjrj, . , . 

— = -AttDM^I - -LL{p + Peff ), (8) 

where D = A A is a positive constant. In the asymptotic regime r — > oo, we have 

M = -An DM 2 (p + p eff ) , (9) 

which clearly demonstrates the increase in mass of the wormhole if p + p e ff < and 
vanishing in the opposite case. 

3 Accretion of viscous phantom energy 

We now study the evolution of mass of WH in two special cases: (a) the constant viscosity; 
and (b) the power law viscosity. 

3.1 Constant bulk viscosity 

For constant viscosity, the evolution of a(t) is given by [10] 

a(t) = a £ % + 7 # £?(t)]*, (10) 



where 

B(t) = exn I - 



B{t) = exp ( 3 -^f] - 1, (11) 



and 7 < 0. The density evolution is given by 



Further, for 7 < the big rip singularity occurs in a finite time at 

-sH 1 -^)- (13) 

Finally, the WH mass evolution is determined by using Eqs. (10)-(12) in (9) to get 



M = M a 



'2£ 0l / \, 2H B{t){i - 1 H ) io& , fa 



l-4nDM ^ln - + - — m + 3 ^ ln V 



-i -1 



(14) 

We can analyze this expression in some asymptotic limits. Assume £ Q is very large i.e. 
£ ^> jH B(t) and £ Q ^> jH while a ~ a a . Hence Eq. (14) reduces to 

M « M [l - 8TTDM H B(t)]-\ (15) 



3 



Using the approximation t ~ £ Q 1 , Eq. (15) becomes 

M w M (l - 137vDM H )-\ (16) 

which finally becomes 

M^M [l + 137rM # + 0(# 2 )]. (17) 

Notice that the current value of Hubble parameter is H a pa 2.3 x 10 _18 s _1 . Due to 
smallness of H Q , its contribution to higher order terms in Eq. (17) is negligible. One can 
see from Eq. (17) that mass of the wormhole increases under the assumption of large 
viscosity. Hence the wormhole is perfectly supported with the viscous phantom energy. 
Similarly, Eq. (14) can be analyzed when £ D <C 1 and £ Q <C \jH B(t)\, hence we get 

T f ( f.. \ / a \ 1 1 ~ ' 

M pa M Q 



In the big rip scenario when a(t) — > oo, the two quantities in curly brackets in Eq. (18) 
will be of the same order of magnitude having opposite sign and hence cancel each other. 
Thus Eq. (18) reduces to 

M pa M . (19) 



3.2 Power law bulk viscosity 



Now we consider the bulk viscosity to possess power law dependence upon density i.e. 
£ = ap s , where a and s are constant parameters. Let us take £ = ap 1 ^ 2 as a special case. 
Then the scale factor evolves as [10] 



a[t) = a a ( 1 



The density of phantom fluid evolves as 

Pit) 



1 - 



3r 2 ( 7 - V3a) 2 V r 
or in terms of critical density p cr as 



-2 



Pit) = Per (l ~ -) 



-2 



The corresponding big rip time r is given by 

2 



T 



-H7 1 . 



(20) 



(21) 



:22) 



(23) 



3(v^3« — 7) ° 

Finally, the mass evolution of WH is determined by using Eqs. (20)- (22) in (9), we get 



M = M Q 



1 - 4itDM < jrp c 



1 - 



+ 3ay/p~ c 



-1 



In a — In a r 



(24) 

In the big rip scenario when t — > r and a — > 00, the two quantities in the square brackets in 
Eq. (24) will be added. This expression gives the growth of the wormhole if the dominant 
energy condition (p + p > 0) is violated. 
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4 Conclusion 

Wormholes are tunnel like topological structures supported by the exotic matter like 
phantom energy. In our model, we have incorporated the viscous pressure along the usual 
isotropic pressure in the accretion model. Our model predicts the growth of wormholes 
by the accretion of bulk viscous phantom energy. For large bulk viscosity, the increase in 
the mass of wormhole is large as compared to small viscosity. 
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